Nonlinear absorption and dispersion in fiber-taper-coupled silicon photonic crystal microresonators 
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A technique is demonstrated which efficiently transfers light between a tapered standard single-mode optical 
fiber and a high-Q, ultra-small mode volume, silicon photonic crystal resonant cavity. Cavity mode quality fac- 
tors of 4.7 x 10 4 are measured, and a total fiber-to-cavity coupling efficiency of 44% is demonstrated. Using this 
efficient cavity input and output channel, the steady-state nonlinear absorption and dispersion of the photonic 
crystal cavity is studied. Optical bistability is observed for fiber input powers as low as 250 /jW, corresponding 
to a dropped power of 100 /jW and 3 fJ of stored cavity energy. A high-density effective free-carrier lifetime for 
these silicon photonic crystal resonators of ~ 0.5 ns is also estimated from power dependent loss and dispersion 
measurements. 

PACS numbers: 



I. INTRODUCTION 

Recently it has been demonstrated that resonant microcavi- 
ties formed in planar photonic crystals (PC) are capable of not 
only confining light to ultra-small optical mode volumes, but 
also of realizing photon cavity lifetimes large enough | jj|2l>|3| 
to enable, for instance, reaching strong-coupling with atomic 
Cs 1 4] or semiconductor quantum dots |5l |a 01 • The en- 
hancement to the local energy density enabled by PC cavi- 
ties is also of significant interest in nonlinear optics, as it re- 
duces the input power required for nonlinear effects, such as 
optical bistability |8, 9]. However, in addition to requiring 
large field enhancements and photon lifetimes, many of the 
proposed applications of PC cavities in quantum, nonlinear, 
and integrated optics, demand the ability to efficiently inter- 
face the PC cavity with external optics. For example, using 
high-Q PC cavities for chip based cavity-QED (cQED) |4] 
or in single -photon sources ilfll . where photon collection is 
an important measure of device performance fTllfl2Hl3ll . re- 
quires an efficient coupling scheme to the sub-micron cavity 
mode. Similarly, proposed applications of microresonators in 
nonlinear-optical switching circuits 1 14] puts a premium on 
reduced input power and efficient optical coupling techniques. 

The difficulty in optically accessing PC cavities is largely 
a result of their ultra-small mode volume and external radia- 
tion pattern, which unlike micropost 1 15] and Fabry-Perot lll6Tl 
cavities, is not inherently suited to coupling with conventional 
free-space or fiber optics. Only recently have low-loss fiber 
coupling techniques to planar photonic crystal waveguides 
been demonstrated experimentally, using on-chip spot size 
converters (~ 1-4 dB loss/port fj^fill'). o ut-of plane diffrac- 
tion gratings (~ 8-10 dB loss/port IU91 I201). and narrow-band 
evanescent coupling (~ 0.1 dB loss/port t2W ). In this paper 
we employ the latter technique to efficiently source and collect 
light from a silicon (Si) high-Q PC cavity via a photonic crys- 
tal waveguide (PCWG) 1 22] which is evanescently coupled to 
a fiber taper I23ll . Using a mode-matched PCWG-PC cavity 
design 1 24] to minimize parasitic loading of the PC cavity, 
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cavities loaded to 60% of critical and maintaining a 2-factor 
close to 4 x 10 4 are demonstrated, and a total fiber-to-cavity 
coupling efficiency of 44% is measured. The utility of this ef- 
ficient fiber coupling technique is then demonstrated in stud- 
ies of the power dependent nonlinear response of the PC cav- 
ity. We observe optical bistability for 100 /jW dropped cavity 
power, and predict sub-nanosecond free carrier lifetimes in the 
cavity. 

An outline of the paper is as follows. We begin in Section 
|n]with a description of the fiber-to-cavity coupling scheme, 
with specific emphasis on issues associated with the complex 
modal properties of a photonic crystal microcavity. In Sec- 
tion |lll| we augment this linear theory by incorporating non- 
linear processes which depend upon the magnitude of the in- 
ternally stored cavity energy. Low power measurements of 
a fabricated silicon photonic crystal waveguide-cavity system 
are presented in Section IIVI where optical fiber taper prob- 
ing is used to both confirm the localized nature of a PC cavity 
resonance and to study the losses in the fiber-cavity system. 
Higher power measurements in which nonlinear effects be- 
come apparent are studied in Section [V] and the model pre- 
sented in Section[ni]is used to estimate the scale of the differ- 
ent nonlinear processes within the silicon PC cavity. Finally, 
a summary is given in Section lvTI 



II. COUPLING SCHEME - THEORY 

An illustration of the coupling scheme is shown in Figs. 
[Qa) and^b). In this scheme, evanescent coupling between an 
optical fiber taper and a PCWG is used to interface with the 
PC chip. Once on the chip, light is guided to a PC microcavity 
at the terminus of the PCWG. Light that is reflected from the 
PC cavity is then recollected into the backward propagating 
fiber taper mode, where it is separated from the forward prop- 
agating input signal using a fiber splitter. In this way a single 
optical fiber is used to both source and collect light from the 
PC microcavity. 

The PC waveguide and PC microcavity studied here are 
formed from a two-dimensional photonic crystal consisting 
of a square lattice array of air holes in an optically thin slab 
waveguide. This PC cavity-waveguide system was previously 
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FIG. 1: (a) Schematic of the fiber taper to PC cavity coupling 
scheme. The blue arrow represents the input light, some of which 
is coupled contradirectionally into the PCWG. The green arrow rep- 
resents the light reflected by the PC cavity and recollected in the 
backwards propagating fiber mode. The red colored region repre- 
sents the cavity mode and its radiation pattern, (b) Illustration of 
the fiber-PC cavity coupling process. The dashed line represents the 
"local" band-edge frequency of the photonic crystal along the waveg- 
uide axis. The step discontinuity in the bandedge at the PCWG - PC 
cavity interface is due to a jump in the longitudinal (£) lattice constant 
(see Fig.|2j. The parabolic "potential" is a result of the longitudinal 
grade in hole radius of the PC cavity. The bandwidth of the waveg- 
uide is represented by the gray shaded area. Coupling between the 
cavity mode of interest (frequency a>o) and the mode matched PCWG 
mode (cowg = <*•()) ls represented by "f Q , coupling to radiating PCWG 
modes is represented by Yj >0 > an d intrinsic cavity loss is represented 
by y. (c,d) Magnetic field profile, calculated using FDTD, of the 
high-Q PC cavity mode and the fundamental TEi PCWG mode, 
respectively. 



studied theoretically in Ref. [24], where the fundamental 
(TEi) mode of the PCWG was designed to mode-match with 
the fundamental (A®) cavity mode. The mode pattern of the 
TEi PCWG mode and the A® cavity mode, calculated using 
the finite-difference time-domain (FDTD) method, are shown 
in Fig. ^c) and Fig. \J^d), respectively. The mode-matched 
cavity acts as a mirror with high modal reflectivity, r (co), ex- 
cept at frequencies of the localized cavity states, where light 
can resonantly tunnel between the PCWG and the cavity. Pre- 
vious studies have shown that the optical fiber-PCWG evanes- 
cent coupler has near unity coupling efficiency (97%) over a 
bandwidth of roughly 15 nm |21], thus providing the neces- 
sary efficient fiber-chip optical interface. Also, in an indepen- 
dent study, the fundamental A® mode of the graded square lat- 
tice PC cavity was measured to have a Q-factor of 4 x 10 4 and 
a mode localization consistent with an effective mode volume 
ofV eff = 0.9(V") 3 lH- 



Here we combine these two complementary components, 
and study the efficiency with which the PC cavity is loaded 
by the mode-matched, fiber-coupled PCWG. We begin with a 
review of some of the key parameters describing the loading of 
a general resonant structure, pointing out specific implications 
for planar PC microcavities. 



As proposed in Ref. |25] in the context of microsphere res- 
onators, the interaction between a PC cavity and an external 
PC waveguide can be described by two key parameters, the 
coupling parameter K and the ideality factor I: 



K = 



(i) 

(2) 



where the cavity mode is characterized by its resonance fre- 
quency co„, its intrinsic photon loss rate in absence of the ex- 
ternal PCWG (Y), and its coupling rates to the fundamental 
TEi mode and higher order (including radiating) modes 
of the external PCWG (Yj>o)- 1 describes the degree of "good" 
loading, via the PCWG TEi mode in this case, relative to the 
total loading of the resonator. K, on the other hand, is the ratio 
of "good" loading to the parasitic and intrinsic loss channels 
of the resonator. 

The coupling parameter K determines the on-resonance 
fraction of optical power reflected by the cavity back into the 
PCWG mode, 



Ro((0 o ) 



(3) 



The remaining fractional power, 1 —R ((0 o ), is absorbed in- 
side the PC cavity or radiated into the parasitic output chan- 
nels. The reflection resonance full-width at half-maximum 
(FWHM) linewidth is given by the sum of the loss rates for 
all of the loss channels of the cavity, 8(0 = "f + E; Yj- From 
R ((a ) and 8co the quality factor of the PC cavity mode due to 
intrinsic and parasitic loss (i.e., those loss channels other than 
the "good" PCWG TEi channel) can be determined, 



Qi+p = 2Q T 



1 



l±VX(co„) 



Qt(1+K), 



(4) 



where the total loaded quality factor is Qj = co„/8co, and 
where the ± corresponds to the under- and over- coupled 
(K Sj 1) loading condition. On resonance, full power trans- 
fer (critical coupling) from the "good" loading channel to the 
resonant PC cavity mode occurs when K = 1 . 

Whereas K determines the amount of power dropped by the 
resonator, the role of / is more subtle. In the case of an internal 
emitter, the collection efficiency (r|o) of emitted photons into 
the "good" loading channel is given by, 
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(5) 



which depends only upon the coupling parameter K. How- 
ever, the cost of obtaining a large collection efficiency is mea- 
sured by the drop in loaded quality factor of the resonant cav- 
ity mode, which can be written in terms of K, I, and Qj as 
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Qi 



K 



m+K)) 



1-— 



no 
/ 



(6) 



3 



Thus, for a given collection efficiency, to maintain a long cav- 
ity photon lifetime, I should be maximized. 

Utilizing a cavity loading method with / ~ 1 is also impor- 
tant for cavity based nonlinear optics. A simple argument can 
be made by studying the stored energy inside a resonant cavity 
for a given input power. One can write for the on-resonance 
internal stored energy U, 



U = (l-R o ((0 o ))^Pi 



AK I-K{\-I)Qi 



co„ 



(1+*) 



Pi (7) 



where f, is the input power in the "good" loading channel. 
The maximum stored energy in the resonator occurs at K mta = 
7/(2-7), giving a peak stored energy f/ max = l(Qi/<a )P u 
which scales directly with I. 

The integrated PC cavity-waveguide design employed here 
has two important features which serve to maximize I: (i) the 
waveguide and cavity modes have similar transverse field pro- 
files (see Figs. Qc-d)) which allows the cavity to be efficiently 
loaded end-on, and (ii) the end-fire PCWG-cavity geometry 
restricts the cavity to a single dominant output channel, in con- 
trast to side-coupled |2], in-line lE3l . and direct taper coupled 
[3] geometries, in which the cavity radiates equally into back- 
ward and forward propagating waveguide modes (bounding 
K < 1 and / < 0.5). In comparison to other microcavity sys- 
tems, the geometry of the PC cavity-waveguide studied here is 
analogous to a Fabry-Perot cavity with a high reflectivity back 
mirror and a lower reflectivity front mirror through which a 
mode-matched input beam sources the cavity. The geome- 
try is also similar to a side-coupled traveling wave resonator, 
such as a microsphere or microdisk that supports whispering- 
gallery type modes, in which the traveling wave resonance ra- 
diates selectively into a single phase matched output channel. 



III. INFLUENCE OF NONLINEAR ABSORPTION AND 
DISPERSION ON CAVITY RESPONSE - THEORY 



Owing to the ultra-small mode volume and long resonant 
photon lifetimes of the PC cavities studied here, the resulting 
stored electromagnetic energy density can be extremely large 
even for modest input powers (< mW), resulting in highly 
nonlinear behavior of the resonant cavity system. In order 
to account for nonlinear effects one may modify the (lin- 
ear) analysis of Section Hi Al bv allowing the various cavity 
and coupling parameters to depend upon the stored cavity en- 
ergy, a reasonably easy quantity to estimate from experimental 
measurements. In this section, relevant nonlinear processes 
are explicitly incorporated into the description of the cavity 
response through use of carefully defined effective modal vol- 
umes and confinement factors appropriate to nonlinear pro- 
cesses in high-index contrast photonic crystal structures. We 
begin with a description of nonlinear absorption, which tends 
to drive the steady-state nonlinear response of the PC cavities 
studied here. 



A. Nonlinear absorption 

Nonlinear absorption adds power dependent loss channels 
to the photonic crystal cavity, degrading the quality factor as 
the internal cavity energy is increased, which in turn modifies 
the coupling efficiency from the PCWG loading channel. This 
effect is incorporated into the formalism presented in Section 
111 Al bv writing the intrinsic cavity loss rate, "f, explicitly in 
terms of its various linear and nonlinear components: 



i{U)= Yrad + Ylin + YlPA (U) + YfCa(^) • 



(8) 



At low power, the "cold cavity" loss rate is given by Yrad an d 
Yiin, which represent loss due to radiation and linear material 
absorption, respectively. Power dependent nonlinear loss is 
given here by Ytpa afl d YfcA' which represent two-photon and 
free-carrier absorption, respectively; other nonlinear absorp- 
tion processes can be included analogously. The coupling pa- 
rameter K and the quality factor Qi+p depend on Y, requir- 
ing the solution of a system of self-consistent equations for U 
in order to determine the on-resonance cavity response for a 
given PCWG input power Pf. 



U 



K(U) 



co„ 



4K(U) Q i+P {U) 
(l+K(U)) 2 co„ 

fo 



(9) 
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Before eqs. (191-11 1> can be solved, explicit expressions for the 
energy-dependent contributions to Y are required. Beginning 
with relations for nonlinear absorption in bulk media, and tak- 
ing into account the complicated geometry of the PC cavity, 
we now derive expressions for y TPA and Yfca- These expres- 
sions can be written in terms of the internal cavity energy, 
known material parameters, and modal parameters which ac- 
count for the mode shape and localization of the PC cavity 
field. 



1. Two-Photon Absorption 

The 1500 nm operating band of the devices studied in this 
work lies in the bandgap of the host silicon material. For the 
doping densities of the p-type silicon membrane used to form 
the PC cavity (p ~ 1 — 3 Q. cm, Na < 10 16 cm -3 ), free-carrier 
absorption due to ionized dopants is small (0C/ e ~ 10~ 2 cm -1 ). 
Two-photon absorption however, is significant! 20. 1271 123. 12911 . 
especially in the highly localized PC cavities. For a given field 
distribution, the (time-averaged) two-photon absorption loss 
rate at position r can be written as 



Yr PA (r) = P'(r)i eo » 2 (r)£ 2 (r), 



(12) 



where E(r) is the amplitude of the complex electric field pat- 
tern E(r) of the resonant mode of the cavity, e e is the per- 
mittivity of free space, and «(r) is the local (unperturbed) re- 
fractive index. The real, physical electric field of the resonant 
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cavity mode can be written in terms of the complex mode pat- 
tern as E(r,f) = (E(r)e- too ' + E* {r)e +ia ° r )/2. The material 
parameter, p', describes the strength of the two-photon ab- 
sorption process, and can be related to the usual two-photon 
absorption coefficient, p, which relates intensity to loss per 
unit length, by p' = (c/n s ,) 2 p, where c is the speed of light 
in vacuum and n g is the group velocity index associated with 
the measurement of p. Typically, for bulk material measure- 
ments where waveguiding is minimal and material dispersion 
is small, n g can be taken to be equal to n. 

In high-index-contrast photonic crystals, E, n, and p' de- 
pend strongly on spatial coordinate r. Equation d!2t describes 
the local two-photon absorption rate; the effective modal two- 
photon absorption rate which characterizes the absorption of 
the entire cavity mode is given by a weighted average of the 
local absorption rate[45] 



Ytpa 



,fYTPA(r)n 2 (r)£ 2 (r)Jr - U 



fn 2 (r)E 2 (r)di 



Vtpa ' 



(13) 



where P' and Vtpa are defined as, 

^_ /P'(r)» 4 (r)E 4 (r)A- 

P JV(r)£ 4 (r)A- 1 ' 

(,/V(r)£ 2 (r)A-) 2 
VjPA - $n\v)E\T)dv ■ (15) 

In a photonic crystal formed by air holes in silicon, p' (r) = p' si 
inside the silicon and p'(r) = in the air, so that eq. Jl 3I > can 
be written as, 



Ytpa - TtpaPs: 



Ttpa 



Vtpa 
^n\v)E\v)dv 
fn*{r)E 4 {r)dr ' 



(16) 
(17) 



for which f si only integrates over the silicon region of the PC 
cavity. 



2. Free-carrier absorption 

Although, as mentioned above, the (linear) free-carrier ab- 
sorption due to the ionized dopants of the silicon layer used for 
the PC cavities in this work is negligible on the scale of other 
losses, two-photon absorption gives rise to a steady-state pop- 
ulation of electron and hole free-carriers far above this equi- 
librium value. Two-photon absorption induced free-carrier ab- 
sorption thus plays a significant role in the silicon PC cavity 
nonlinear response. At position r in the cavity, assuming a 
simple Drude model, the optical loss rate due to free-carrier 
absorption is 



Yfca = a'(r)jV(r) 



(18) 



where a' is related to the material dependent free-carrier 
cross-section, a, by a' = o(c/n g ), and N(r) is the free- 
carrier density. In silicon it has been demonstrated 



experimentally 1 30] that this model correctly describes absorp- 
tion by both electrons and holes, albeit with unique values of 
G' e h for each carrier type. Here we let N represent the number 
of electron-hole pairs|46], and take a' = <3' e + & h . 

In general, the derivation of the free-carrier density for a 
given two-photon absorbed power distribution requires a mi- 
croscopic theory taking into account carrier diffusion, carrier- 
carrier scattering effects (Auger recombination for instance), 
and in the highly porous PC cavities, local surface recombi- 
nation effects. In lieu of such an analysis, we approximate 
the free-carrier density distribution by considering the local 
two-photon absorbed power, 

tpTPA(r) 



N(t) 



2h&„ 



(19) 



where x is a free-carrier lifetime, and /?tpa(i") is the local ab- 
sorbed power density due to two-photon absorption, 



Ptpa(i-) = ie n 2 (r)£ 2 (r)y T p A (r). 



(20) 



Equation ( I19> neglects non-local effects due to spatial carrier 
diffusion by assuming that N(r) depends only on the power 
absorbed at position r; however, it does correlate regions of 
strong two-photon absorption with high free-carrier density. 
Also, since x generally depends on N and on the proximity to 
surfaces, x will have a spatial dependence within the cavity. 
We neglect this effect here, and let x represent an effective 
free-carrier lifetime for all the carriers in the cavity region| 47 ] . 
Combining eqs. (II 81 , ( 1191 and ( I20> . an effective modal free- 
carrier absorption rate can be written as 



Yfca 



^,/(^(r)l £( y(r)£ 2 (r)YTPA(r))« 2 (r)£ 2 (r)^r 
Jn 2 (r)E 2 (r)dr 



(21) 

Substituting eq. (II 2b for Ytpa(f), the modal loss rate due to 
free-carrier absorption in the porous silicon photonic crystals 
considered here can be written as 

^ / to' p' U 2 
Yfca-" 1 



FCA 



V 2fm v 2 



FCA 



(22) 



with effective confinement factor and mode volume defined as 
J si n 6 (r)E 6 (r)dr 



FCA 



Jn 6 (r)E 6 (r)dr 
(Jn 2 (r)E 2 (r)dr 



(23) 



(24) 



FCA fn 6 (r)E 6 (r)dr 

Equations i22i and H6i represent the total loss rate of pho- 
tons from the cavity due to free-carrier and two-photon ab- 
sorption, respectively. These expressions depend only on ma- 
terial parameters, modal confinement factors Ttpa.fca, effec- 
tive mode volumes Vtpa.fca, and the internal cavity energy 
U. The modal parameters take account of the non-trivial ge- 
ometry and field distribution of the cavity mode, and can be 
determined for a given mode from finite difference time do- 
main simulations. Including expressions (I22> and (II 61 in "f, 
eqs. (191-11 1> can be solved iteratively for Qi + p and K, which 
characterize the on-resonance nonlinear response of the cavity 
for a given input power. 
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B. Nonlinear and thermal dispersion 

In addition to modifying the cavity quality factor, large cav- 
ity energy densities also modify the refractive index of the 
cavity, resulting in a power dependent resonance frequency. 
Here we consider the role of the Kerr effect, free-carrier dis- 
persion, and heating due to linear and nonlinear absorption, 
on the dispersive response of the PC cavity. The refractive in- 
dex shift induced through the processes considered here is a 
function of both space and internal cavity energy. The renor- 
malization of the resonant cavity frequency, resulting from 
small local perturbations in the refractive index, can be ap- 
proximated using first order perturbation theory as 



Aco (C7) 



co„ 



= -An(U), 



(25) 



where the normalized modal index shift, An(U), is given by 
an average of the (normalized) local refractive index shift 
An(r)/n(r), 



An(U) 



J(ffi)"W(r)rfr 
Jn 2 (r)E 2 {r)dr ' 



(26) 



This energy dependent frequency shift, together with the en- 
ergy dependent loss described in Section IIII Al modifies the 
Lorentzian frequency dependence of the cavity response: 



/? (to) = l- 



4K{U) 



(5co/2) 2 



(1 +K(U)) 2 (co - co„ - Aco„(t/)) 2 + (Sco(t/)/2) 2 ' 

(27) 

For a given input power P, and frequency CO, U is given by 



U = - — = 1 -R {(0)) — Pi, 

Ji+P co 



(28) 



where Pj = ( 1 — R a (co) )P, is the frequency dependent dropped 
power in the resonant cavity. 

For input powers sufficient to shift Aco„ > v / 3§co/2, the fre- 
quency response described by eq. 1271 is bistable, and can be 
exploited for applications including temperature locking and 
optical switching In order to solve eq. (1271 for 

the cavity response, it is necessary to derive expressions for 
each of the constituents of An as a function of U. We begin 
with the Kerr effect. 



1. Kerr effect 

The time-averaged local index shift induced by the Kerr ef- 
fect is 

A» K en-(r) = ^(r)i £o n 2 (r)£ 2 (r), (29) 

where n' 2 {r) is a material parameter, and is related to the usual 
«2 coefficient relating intensity to refractive index shift 1 34] by 
n' 2 = (c/n s ,)«2- In a silicon PC cavity the normalized modal 
index change due to the Kerr effect can be written as|48] 



An K err(£/) = n 2si 



"Si 



Vkerr J ' 



(30) 



with n' 2 si and «si the Kerr coefficient and linear refractive in- 
dex of Si, respectively. As both the Kerr effect and two photon 
absorption (TPA) share the same dependence on field strength, 
the confinement factor and effective mode volume associated 
with the Kerr effect are equal to those of TPA, 

TKerr = Ttpa (31) 
VKerr = VtPA- (32) 

2. Free-carrier dispersion 

Dispersion due to free-carrier electron-hole pairs is given 
by|49] 



AnpcD(r) = -C(r)iV(r), 



(33) 



where £(r) is a material parameter with units of volume. Fol- 
lowing the derivation of Y F ca> the normalized modal index 
change is 



A«FCD(t/) 



r FCD Z-cCsiPsi u 2 



«Si V 2h ™o V? 



FCD 



with 



(34) 



(35) 
(36) 



Tfcd = Tfca 
Vfcd = Vfca- 

3. Thermal dispersion 



It is also necessary to consider the effect of thermal heating 
due to optical absorption on the refractive index of the PC 
cavity. The normalized modal index shift is given by 



A/ith 



n^{r)6T{T))f?{T)E 2 {T)dv 



J ' n 2 {v)E 2 {v)dr 



(37) 



Here AT(r) is the local temperature change due to the ab- 
sorbed optical power density within the cavity, and dn /dT is 
a material dependent thermo-optical coefficient. Neglecting 
differences in the spatial distributions! 50] of the contributions 
to AT(r) from the various absorption processes, and assuming 
that AT(r) scales linearly with absorbed power density for a 
fixed spatial heating distribution, the modal index shift can be 
written as 



, , r t h / dnvi dT , „ 



(38) 



where, 



P,bs(U) = (Yim +?tpa(^) +1fca(U 2 )) U. (39) 



r t h = J si n rE dxj J n E dv is a confinement factor which ac- 
counts for the fact that only the semiconductor experiences 
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an appreciable index shift, and dT / JP a bs is the thermal resis- 
tance of the PC cavity which relates the mean modal temper- 
ature change to the total absorbed power. In what follows we 
lump these two factors together, yielding an effective thermal 
resistance of the PC cavity. 

From eqs. j30L 1341 . and j38> . the total modal index change 
and corresponding resonance frequency shift can be deter- 
mined as a function of cavity energy. The nonlinear lineshape 
described by eq. \21\ can then be calculated iteratively as a 
function of input power when combined with the power de- 
pendent loss model of sub-section IIIIAI This is used below 
in Section [V] to estimate the scale of the different nonlinear 
processes in silicon PC microcavities. 

As a final comment we note that the above analysis has 
assumed a steady-state optical, carrier, and thermal distribu- 
tion, whereas of significant practical interest for applications 
such high speed switching is the transient response of such 
structures. Although the Kerr nonlinearity, two-photon ab- 
sorption, free-carrier absorption, and free-carrier dispersion 
all depend on the electronic structure of the semiconductor 
material, the sub-micron geometry typical of photonic crys- 
tals can also play an important role. For example, the sur- 
faces introduced by the slab and air hole geometry of planar 
PC cavities can significantly modify the free-carrier lifetime x 
compared to that in bulk material 1 2 ill 1351 13611 . Similarly, the 
thermal response time scales inversely with the spatial scale 
of the optically absorbing region, and depends upon the ge- 
ometry and material dependent thermal properties of a given 
structure 13411 . Although not the focus of the work presented 
here, an inkling of these effects is seen in the sub-nanosecond 
estimated effective free-carrier lifetime in the silicon PC cav- 
ity studied below. 



IV. EFFICIENT COUPLING INTO A HIGH-Q CAVITY 
MODE 

The integrated PC cavity-PCWG devices described in Sec- 
tion|H]were fabricated in an optically thin layer (thickness 340 
nm) of silicon as described in Ref. |37]. A typical device is 
shown in Fig. which also shows the regions in which un- 
patterned silicon was removed to allow taper probing of the 
PC devices. In addition to isolating the PC devices on a mesa 
of height ~ 10 /jm, a trench extending diagonally from the 
cavity was defined. This allowed the cavity to be probed di- 
rectly by the fiber taper as in Ref. |3J]. Indirect, yet much more 
efficient coupling to the PC cavity, was performed using the 
scheme described in Sec.[|1]by aligning the fiber taper along 
the axis of the PCWG and coupling through the PCWG into 
the PC cavity. The fiber taper was mounted in a "u"-shape 
on a DC motor z-stage with 50 nm encoder resolution, pro- 
viding accurate vertical placement of the taper above the Si 
chip. In the measurements described below the taper was held 
fixed in the x-y plane parallel to the Si chip surface, and in- 
plane positioning of the Si chip was performed using a pair of 
DC motor stages with similar encoder resolution. A rotation 
and goniometer stage, mounted to the x-y stages, was used to 
align the in-plane and vertical angle of the sample relative to 



the taper, respectively. The entire sample and taper system 
was enclosed in an acrylic box to reduce taper fluctuations 
due to air currents in the room. Further details of the fiber ta- 
per fabrication and mounting geometry can be found in Refs. 




FIG. 2: SEM image of an integrated PCWG-PC cavity sample. The 
PC cavity and PCWG have lattice constants A ~ 430 nm, A x ~ 430 
nm, and A z ~ 550 nm. The surrounding silicon material has been 
removed to form a diagonal trench and isolated mesa structure to 
enable fiber taper probing. 




Fiber taper orientations: (i) (ii) 

FIG. 3: (a) Illustration of the device and fiber taper orientation for (i) 
efficient PCWG mediated taper probing of the cavity, and (ii) direct 
taper probing of the cavity, (b) Normalized depth of the transmission 
resonance (AT) at X ~ 1589.7, as a function of lateral taper dis- 
placement relative to the center of the PC cavity, during direct taper 
probing (taper in orientation (ii)). 



A fiber-coupled swept wavelength (1565 - 1625 nm) laser 
source was used to measure the wavelength dependent for- 
ward transmission (T) through the fiber taper, and a fiber split- 
ter was used to monitor the signal in the backward propagat- 
ing fiber taper mode (R). With the taper aligned with the 
etched trench, the spectral and spatial properties of the PC 
cavity modes were probed directly (Fig. |3Ja), taper position 
(ii)), as in Ref. |3]. The trench prevents the fiber taper from 
interacting with the unpatterned silicon, and light is coupled 
directly from the fiber taper into the high-Q PC cavity modes. 
Although this coupling is inefficient (AT =1-10%, / <C 1), it 



7 




1584 . 1586 ' 1588 1590 1592 1594 



Wavelength (nm) 




Wavelength (nm) Taper height g (jim) 

FIG. 4: (a) Measured reflected taper signal as a function of input 
wavelength (taper diameter d ~ 1 //m, taper height g = 0.80 /jm). 
The sharp dip at X ~ 1589.7 nm, highlighted in panel (b), corre- 
sponds to coupling to the Aj cavity mode, (c) Maximum reflected 
signal (slightly detuned from the Aj resonance line), and resonance 
reflection contrast as a function of taper height. The dashed line at 
AR = 0.6 shows the PCWG-cavity drop efficiency, which is indepen- 
dent of the fiber taper position for g > 0.8 /jm. 



allowed the frequency of the A" cavity mode to be indepen- 
dently determined. 

In the device studied here, when the taper was aligned with 
the trench and positioned ~ 500 nm above the cavity, a sharp 
dip in T was observed at a wavelength of X ~ 1589.7 nm. 
It was confirmed that this was due to coupling to a localized 
cavity mode by studying the depth of the resonance as the 
taper was displaced laterally (|Ax| > 0) relative to the center of 
the PC-cavity. The measured normalized resonance depth as 
a function of taper displacement is shown in Figure[3Jb), and 
has a halfwidth of 480 nm, consistent with previous studies of 
the localized Aj cavity modey|]. 

The fiber taper was then aligned above and parallel to the 
PCWG (Fig. a), taper position (i)). At taper-PCWG phase- 
matching wavelengths, T decreases resonantly as power is 
coupled from the taper into the PCWG; coupling to the TE] 
PCWG mode was verified by studying the dispersive and spa- 
tial properties of the coupling, as in Ref. 1 37 ] . The fiber taper- 
PCWG coupling bandwidth was adjusted to overlap with the 
wavelength of the A^ cavity mode using two mechanisms. 
Coarse tuning was obtained by adjusting, from sample to sam- 
ple, the nominal hole size and longitudinal lattice constant 
(A z ) of the PCWG. Fine tuning of the coupler's center wave- 
length over a 100 nm wavelength range was obtained by ad- 
justing the position, and hence diameter, of the fiber taper re- 
gion coupled to the PCWG l37ll . Different degrees of cavity 
loading were also studied by adjusting the number of periods 
(9-11) of air holes between the center of the PC cavity and 
the end of the PCWG. In the device studied below (shown in 
Fig. |2J the PC cavity was fabricated with 9 periods on the side 
adjacent the PCWG and 18 periods on the side opposite the 
PCWG . 

Figure Ua) shows the normalized reflected fiber signal, R, 



for a taper diameter d ~ 1 /jm, which aligns the taper-PCWG 
coupler bandwidth with the Aj PC cavity mode wavelength. 
This signal is normalized to the taper transmission in absence 
of the PCWG, and since light passes through the taper-PCWG 
coupler twice, is given by R — r\f vg R , where T[ wg is the taper- 
PCWG coupling efficiency. Note that both R and r\ wg are 
frequency dependent. In Figure Ufa), the peak in R around 
X ~ 1590 nm corresponds to the phase-matched point of the 
fiber taper and the TEi PCWG mode. From the peak value of 
R max = 0.53, a maximum taper-PCWG coupling efficiency of 
Twg ~ 73% is estimated, where the off-resonant R is taken to 
be unity. This value is lower than the 97% obtained in pre- 
vious work 1211 due to coupling to additional higher-order 
(normal to the Si slab) PCWG modes which interfere with 
the coupling to the fundamental TEi PCWG mode for strong 
taper-PCWG coupling. This can be avoided in future devices 
by increasing the nominal PCWG hole size relative to that in 
the PC-cavity or reducing the Si slab thickness, effectively 
freezing out the higher-order PCWG modes l37Tl . 

The sharp dip in reflection at X ~ 1589.7 nm, shown in de- 
tail in Fig. Elb), corresponds to resonant excitation of the A^ 
PC cavity mode, as confirmed by the direct fiber probing of 
the cavity described above. The other broad features in R cor- 
respond to weak Fabry-Perot effects of the PCWG, and weak 
interference between the TE- 1 mode and higher order PCWG 
modes. The reflected fiber taper signal as a function of taper- 
PCWG gap height, g, is shown in Fig. Etc). For g > 0.8 
/jm, R mii x increases with decreasing g as the coupling from the 
fiber taper to the TEi PCWG mode becomes stronger. The re- 
flection contrast, AR = 1 -R ((d ) = (R max —R((Oo))/RrBax, 
remains constant, since the PCWG-cavity interaction is inde- 
pendent of the fiber taper to PCWG coupling. For smaller 
taper-PCWG gap heights, g < 0.8 /mi, fiber taper coupling 
into higher order PCWG modes and radiation modes becomes 
appreciable, and R max decreases for decreased taper height. 
The corresponding increase in AR seen in Fig. |4jc) is a re- 
sult of interference between the TEi mode and higher-order 
PCWG modes which are excited and collected by the taper, 
and is not a manifestation of improved coupling between the 
TEi PCWG mode and the A^ PC cavity mode. Direct coupling 
between the taper and the cavity is negligible here. 

From a Lorentzian fit to the Aj cavity resonance dip in 
R o {03), the normalized on-resonance reflected power is esti- 
mated to be R o (03 o ) = 0.40, corresponding to an undercou- 
pled K — 0.225. The loaded quality factor as measured by 
the reflected signal linewidth is Qt = 3.8 x 10 4 . Substituting 
these values into eq. @ gives the cavity mode quality factor 
due to parasitic loading and intrinsic losses, Qi+p = 4.7 x 10 4 . 
Previous measurements of similar PC cavity devices without 
an external PCWG load yielded intrinsic quality factors of 
4 x 10 4 yl, strongly indicating that the parasitic loading of 
the PC cavity by the PCWG is minimal, and / ~ 1 for this PC 
cavity-waveguide system. The high ideality of this coupling 
scheme should be contrasted with previous direct taper mea- 
surements of the PC cavities |3], whose coupling was limited 
to a maximum value K = 0.018, with an ideality of / ~ 0.035 
(corresponding to a resonance depth of 7%, Qt = 2.2 x 10 4 ). 

The efficiency of power transfer from the fiber taper into 
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the PC cavity is given by T|,„ = T\ H , g AR w 44%. This corre- 
sponds to the total percentage of photons input to the fiber 
taper which are dropped by the PC cavity. Based upon these 
measurements, in the case of an internal cavity source such 
as an atom or a quantum dot, the efficiency of light collec- 
tion into the fiber taper for this PC cavity system would be 
TW = T| wg r|o w 13% (r|o m 18%). As the tapers themselves 
are of comparably very low loss, with typical losses associ- 
ated with the tapering process less than 10%, these values ac- 
curately estimate the overall optical fiber coupling efficiency. 
Finally, note that previous measurements of near-ideal cou- 
pling between the fiber taper and PCWG 12111 indicate that 
by adjusting the PCWG as described above to improve X\ wg , 
T[i n and r\ out can be increased to 58% and 18%, respectively. 
More substantially, adjustments in the coupling parameter K 
towards over-coupling by decreasing the number of air-hole 
periods between the PC cavity and the PCWG can result in 
significant increases in T|,„ and T[ out with minimal penalty in 
loaded Q-factor for / ~ 1. 



V. NONLINEAR MEASUREMENTS 

The nonlinear response of the PC cavity was studied by 
measuring the dependence of the reflected signal lineshape on 
the power input to the PCWG. Figure [5] shows wavelength 
scans of the cavity response R a for varying power, P r , input 
to the fiber taper. Each scan was obtained by dividing the 
normalized reflected signal, R, by the slowly varying taper- 
PCWG coupler lineshape, T|jL ((*)). In all of the measure- 
ments, the fiber taper was aligned near the optimal taper- 
PCWG coupling position, and the wavelength of the laser 
source was scanned in the direction of increasing X. P t was 
determined by taking taper insertion loss into account, and 
measuring the taper input power with a calibrated power me- 
ter. 

Increasing the power in the fiber taper, and consequently the 
PCWG, results in three readily observable changes in R : (i) a 
decrease in the resonance contrast, AR ( „ (ii) a shift Aco„ in the 
resonance frequency C0 o , and (iii) broadening and asymmet- 
ric distortion of the resonance lineshape, eventually leading to 
a "snap" in the reflection response characteristic of bistabil- 
ity IB8I1 . Here we use the theory presented in Section ITTT1 to 
show that these features are due to nonlinear absorption and 
dispersion in the PC cavity. 

Figure^a) shows P ( /, the on-resonance power dropped into 
the PC cavity, as a function of p, the power incident on the 
cavity from the PCWG. P c / is measured from P c { = AP (p)p, 
and Pi is related to the taper input power by P, = r[ wg (G} )P t . 
For small p, P/ increases with a constant slope equal to the 
"cold cavity" value of AR U = 0.60 measured in Section Hvl 
For larger P„ P c / becomes sub-linear versus P,- as loss due 
to nonlinear absorption becomes appreciable compared to the 
other loss channels of the PC cavity. In the context of the anal- 
ysis of Section||ll] Y increases with increasing P t /, degrading 
K, and decreasing AP„ (for K < 1). From the "cold cavity" 
T|o and Qt measured in the previous section, the power de- 
pendent Qi+p(Pj) can be extracted from AP„(p) through the 
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FIG. 5: (a) Measured cavity response as a function of input wave- 
length, for varying PCWG power (taper diameter d ~ 1 /jm, taper 
height g = 0.80 /mi). 
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FIG. 6: (a) Power dropped (P£) into the cavity as a function of power 
in the PCWG (Pi). The dashed line shows the expected result in 
absence of nonlinear cavity loss, (b) Resonance wavelength shift as 
a function of internal cavity energy. Solid blue lines in both figures 
show simulated results. 



relation: 

Q l+P (Pi) = KiARoiPt)) ^^ . (40) 

Equation d40l is useful for powers where nonlinear effects dis- 
tort the Lorentzian lineshape, and X /5X is not an accurate 
measure of Qr(Pi)- Using eqs. J40i and Q, the internal cav- 
ity energy, U, can be calculated from P, and AP G . 

Figure|5Jb) shows a plot of the measured l5lll AX , the res- 
onance wavelength shift, as a function of U. This plot has 
several noteworthy properties. First, the wavelength shift is 
nonlinear in U, indicating that nonlinear processes such as 
free-carrier dispersion and heating through nonlinear absorp- 
tion must be taking effect. Also, for small U the resonance 
wavelength is seen to blue shift. In the 1550 nm wavelength 
band of operation both drisi/dT and «2,Si are > 0, while 
d(Anpco)/dU < 0, indicating that free-carrier dispersion is 
the dominant dispersive process at low input powers. For 
U > 0.34 fJ (P c i > 10 /jW), the resonance wavelength begins 
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to red shift, indicating that thermal or Kerr effects dominate 
for large internal cavity energy. Also, note that for a stored 
cavity energy as low as U ~ 3 O (Pd ^100 /jW) the cavity 
response is bistable with AA,„ = 35 pm ~ \/38A,/2. 

In order to estimate the contributions of the various nonlin- 
ear processes to the effects discussed above, the absorptive, 
P c i{Pi), and dispersive, AA,(P,), data were fit using the model 
presented in SectionlliTl Specifically, eqs. J9I1 1 1 were solved 
for Pd and U as a function of P,, and eq. (I25i was used to 
calculate Ak . The free parameters in this model were taken 
as: (i) the effective free-carrier lifetime, x, (ii) the effective 
thermal resistance of the PC cavity, T^dT /c/P a bs, and (iii) the 
fraction of the "cold cavity" loss which is due to linear ab- 
sorption (as opposed to radiation), r|ij n = Yii n / (Ylin +Yrad)- The 
material and modal constants used are listed in Table|I] 

As has been observed in studies of silicon optical 
waveguides 1 28], we find that a strong dependence of x on car- 
rier density is required for our model to accurately reproduce 
both the dispersive and absorptive data represented in Figures 
0a) and (b). In order to account for a carrier density depen- 
dent lifetime in our model the following procedure was used. 
With FfadT / dP a b s and T|r m held fixed, x(P,) was determined 
for each input power from a least squares fit to AX (Pi) and 
Pd{Pi)\52]. This procedure was repeated for a range of values 
for TfadT I dP a b s and r|ij n . For a fixed value of T|ii n , the fits 
were robust in T^iir / dP&s with the sum of the least square 
residual of x(f,) clearly minimized for an optimal value of 
TfadT I dPabs- This procedure, however, was only found to 
constrain T|ii n > 0.15. Within this range of r\n n the quality of 
the fits does not change significantly, with the optimal func- 
tional form of x changing slightly and the optimal value of 
r±dT /dPabs varying between ~ 15 — 35 K/mW. Based on es- 
timates of T|ii n from studies of loss in silicon microdisk res- 
onators fabricated using the same SOI wafers and the same 
processing techniques 1 39], and by comparing the etched sur- 
face area seen by the PC cavity mode to that seen by a mi- 
crodisk mode, we chose to use r\\ m ~ 0.40 for the PC cavity. 
With this value of r|ii n the optimal value of the effective cavity 
thermal resistance, T^dT /<fP a bs, was found to be 27 K/mW, 
of the same order of magnitude as the result calculated in Ref. 
ll40ll for a similar membrane structure. Finally, the point-by- 
point least-squared optimum values of x(P,) were then fit with 
a smooth curve of functional form t - 1 = A + BN a as a func- 
tion of the effective free-carrier density N. Using this x(A^) 
fit, smooth fits to measured Pd{Pi) and AA,(P,) were obtained, 
shown as solid blue lines in Fig. [6] 

The various components of the total cavity loss rate and 
resonance shift based upon the above fits to the measured data 
are shown in Fig. It can be seen that although TPA does not 
dominate the PC cavity response, the free-carriers it generates 
and the resulting free carrier dispersion and absorption drive 
the nonlinear behavior of the silicon PC cavity at low and high 
input powers, respectively. The fit effective free-carrier life- 
time, shown in Fig. [8] shows similar characteristics to that 
obtained by Tsang, et al. l28ll . demonstrating a significant 
fall-off in x for large N, but with a smaller saturated lifetime. 
Both the pronounced decay in x and the low ~ 0.5 ns value of 



TABLE I: Fixed parameters used in the model. 
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Value 




Units 


Source 


Vtpa 


4.90 




(V«s,) 3 


FDTD(b) 


Vfca 


3.56 




(V«s.) 3 


FDTD(b) 


Ftpa 


0.982 






FDTD(b) 


Ffca 


0.997 






FDTD(b) 


"Si 


3.45 






[22,41,] 




14.5 x 10" 


-22 


m 2 




& 


8.8 x 10- 


28 


m 3 


[20,44] 




4.6 x 10- 


28 


m 3 


[2&41J 


»2,Si 


4.4 x 10" 


18 


nvUw- 1 


[27J 


Psi 


8.4 x 10- 


12 


m-W- 1 


[27J(a) 


dnsi/dT 


1.86 x 10 


-4 


K- 1 


[42] 



(a) Average of the two quoted values for Si(l 10) and Si (111). 

(b) Calculated from FDTD generated fields of the Aj cavity mode 
of the graded square lattice cavity studied here. 
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FIG. 7: (a) Simulated effective quality factors for the different PC 
cavity loss channels as a function of power dropped into the cavity, 
(b) Contributions from the modeled dispersive processes to the PC 
cavity resonance wavelength shift as a function of power dropped 
into the cavity. (Simulation parameters: r|ii n ~ 0.40, T t hdT / dP. d bs = 
27 K/mW, t- 1 ~ 0.0067 + (1.4 x 10- 7 )7Y°- 94 where N has units of 
cm -3 and x has units of ns). 



the high-carrier density free-carrier lifetime are significantly 
different from that found in bulk Si, and are most likely re- 
lated to carrier diffusion and surface effects owing to the ex- 
tremely large surface-to-volume ratio of the PC cavity, the 
small length scales involved (~ 200 nm feature size), and the 
small size scale of the optical mode|35]. This small effective 
free-carrier lifetime is consistent with other recent experimen- 
tal results of highly porous silicon optical structures l3ll l36ll . 
It should, however, be noted that the bulk Si TPA coefficient 
was used in modelling the nonlinear response of the PC cav- 
ity, which given the above comments may not be accurate due 
to surface modification of TPA. As the effects of free-carrier 
lifetime and two-photon absorption on the behaviour of the 
dispersive and absorptive nonlinear response of the PC cavity 
are somewhat intertwined, further studies will be necessary to 
concretely separate these two phenomena in porous Si struc- 
tures such as the photonic crystals of this work. 
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FIG. 8: Dependance of free-carrier lifetime on free-carrier density 
(red dots) as found by fitting &k (Pj) and Pd{Pi) with the constant 
material and modal parameter values of Table and for effective 
PC cavity thermal resistance of T t hdT /dP a bs = 27 K/mW and linear 
absorption fraction T|ij n = 0.40. The solid blue line corresponds to a 
smooth curve fit to the point-by-point least-squared fit data given by 
T" 1 ~ 0.0067 + (1.4 x 10- 7 )N°- 94 
1 is in ns. 



where N is in units of cm 3 and 



Q ultra-small mode volume PC cavities. By employing a 
PCWG which supports a mode which is simultaneously spa- 
tially mode matched with the high-Q PC cavity mode of in- 
terest and phase matched with an optical fiber taper, effi- 
cient fiber-PCWG-cavity coupling is enabled. A total fiber- 
to-cavity coupling efficiency of 44% is demonstrated, which 
for the case of an internal cavity emitter corresponds to a ra- 
diated photon collection efficiency of 13%. These values are 
not fundamental limits of this technique, and can be improved 
through fine tuning of the photonic crystal. The efficiency 
of this fiber coupling method was then exploited to probe the 
steady-state nonlinear optical properties of the PC cavity. The 
effect of two-photon absorption, free-carrier absorption and 
dispersion, Kerr self-phase modulation, and thermo-optic dis- 
persion, on the response of the PC cavity was considered. Op- 
tical bistability at fiber input powers of 250 /jW was observed, 
and a free-carrier lifetime with high carrier density value as 
low as ~ 0.5 ns is inferred from nonlinear absorptive and dis- 
persive measurements of the PC cavity. Along with applica- 
tions to nonlinear optics, this optical fiber evanescent coupling 
based approach should be useful for future experiments in in- 
tegrated micro-optics with photonic crystals, and in particular 
to cQED systems employing photonic crystals, where quan- 
tum computing and communication protocols demand high 
optical fidelity. 



VI. SUMMARY 

In conclusion, we have demonstrated an optical fiber cou- 
pling scheme to efficiently source and collect light from high- 



The authors thank Matthew Borselli for helpful discussion. 
This work was supported by the Charles Lee Powell founda- 
tion. KS would also like to thank the Hertz foundation for its 
support through a graduate fellowship. 



[1] K. Srinivasan, P. Barclay, O. Painter, J. Chen, A. X. Cho, and 
C. Gmachl, Appl. Phys. Lett. 83(10), 1915 (2003). 

[2] Y. Akahane, T. Asano, B.-S. Song, and S. Noda, Nature 425, 
944 (2003). 

[3] K. Srinivasan, P. E. Barclay, M. Borselli, and O. Painter, Phys. 

Rev. B 70, 08 1 306(R) (2004). 
[4] B. Lev, K. Srinivasan, P. E. Barclay, O. Painter, and H. Mabuchi, 

Nanotechnology 15, S556 (2004). 
[5] T. Yoshie, A. Scherer, J. Hendrickson, G. Khitrova, H. Gibbs, 

G. Rupper, C. Ell, O. Shchekin, and D. Deppe, Nature (London) 

432, 200 (2004). 

[6] J. Reithmaier, G. Sek, A. Loffler, C. Hofmann, S. Kuhn, S. Re- 
itzenstein, L. Keldysh, V. Kulakovskii, and T. Reinecke, Nature 
(London) 432, 197 (2004). 

[7] E. Peter, P. Senellart, D. Martrou, A. Lemaitre, and J. Bloch, 
quant-ph/04 11076 (2004). 

[8] A. R. Cowan and J. F. Young, Phys. Rev. E 68, 046606 (2003). 

[9] M. Soljacic, S. G. Johnson, Y. Fink, and J. D. Joannopoulos, 
Phys. Rev. E 66(5), 055601 (2002). 
[10] J.-M. Gerard, Solid State Cavity-Quantum Electrodynamics 
with Self-Assembled Quantum Dots (Springer- Verlag, Ger- 
many, 2003), pp. 269-314. 
[11] G. Brassard, N. Lutkenhaus, T. Mor, and B. Sanders, Phys. Rev. 
Lett. 85(6), 1330 (2000). 



[12] E. Knill, R. Lafiamme, and G. Millburn, Nature (London) 

409(6816), 46 (2001). 
[13] P. Kok, H. Lee, and J. P. Dowling, Phys. Rev. A 66, 063814 

(2002). 

[14] M. F. Yanik, S. Fan, and M. Soljacic, Appl. Phys. Lett. 83(14), 
2739 (2003). 

[15] M. Pelton, C. Santori, J. Vuckovic, B. Zhang, G. Solomon, 
J. Plant, and Y. Yamamoto, Phys. Rev. Lett. 89, 299602 (2002). 

[16] J. McKeever, A. Boca, A. D. Boozer, R. Miller, J. R. Buck, 
A. Kuzmich, and H. J. Kimble, Science 303(5666), 1992 
(2004). 

[17] S. J. McNab, N. Moll, and Y. A. Vlasov, Opt. Expr. 11(22), 
2927 (2003). 

[18] M. Notomi, A. Shinya, S. Mitsugi, and H.-Y. Ryu, Opt. Expr. 

12(8), 1551 (2004). 
[19] W. Bogaerts, D. Taillaert, B. Luyssaert, P. Dumon, J. V. Camp- 

enhout, P. Bienstman, R. Baets, V. Wiaux, and S. Beckx, Opt. 

Expr. 12(8), 1583 (2004). 
[20] A. R. Cowan, G. W. Rieger, and J. F. Young, Opt. Expr. 12(8), 

1611 (2004). 

[21] P. E. Barclay, K. Srinivasan, M. Borselli, and O. Painter, Opt. 

Lett. 29(7), 697 (2004). 
[22] M. Loncar, D. Nedeljkovic, T. Doll, J. Vuckovic, A. Scherer, 

and T. P. Pearsall, Appl. Phys. Lett. 77(13), 1937 (2000). 



11 



[23] J. Knight, G. Cheung, F. Jacques, and T. Birks, Opt. Lett. 

22(15), 1129 (1997). 
[24] P. E. Barclay, K. Srinivasan, and O. Painter, J. Opt. Soc. Am. B 

20(11), 2274 (2003). 
[25] S. M. Spillane, T. J. Kippenberg, O. J. Painter, and K. J. Vahala, 

Phys. Rev. Lett. 91(4), 043902 (2003). 
[26] S. Y. Lin, E. Chow, S. G. Johnson, and J. D. Joannopoulos, Opt. 

Lett. 26(23), 1903 (2001). 
[27] M. Dinu, F. Quochi, and H. Garcia, Appl. Phys. Lett. 82(18), 

2954 (2003). 

[28] T. Liang and H. Tsang, Appl. Phys. Lett. 84(15), 2745 (2004). 
[29] K. Kanamoto, S. Lan, N. Ikeda, Y. Tanaka, Y. Sugimoto, 

K. Asakawa, and H. Ishikawa, IEICE Trans. Electron. E87- 

C(7), 1142 (2004). 
[30] R. A. Soref and B. R. Bennett, IEEE J. Quan. Elec. 23(1), 123 

(1987). 

[31] V. R. Almeida, C. A. Barrios, R. R. Panepucci, and M. Lipson, 

Nature (London) 431(7012), 1081 (2004). 
[32] V. R. Almeida and M. Lipson, Opt. Lett. 29(20), 2387 (2004). 
[33] T. Carmon, L. Yang, and K. J. Vahala, Opt. Expr. 12(20), 4742 

(2004). 

[34] R. W. Boyd, Nonlinear Optics (Academic Press, San Diego, 

CA, 2003), 2nd ed. 
[35] R. Claps, V. Raghunathan, D. Dimitropoulos, and B. Jalali, Opt. 

Expr. 12(12), 2774 (2004). 
[36] H. Rong, A. Liu, R. Nicolaescu, M. Paniccia, O. Cohen, and 

D. Hak, Appl. Phys. Lett. 85(12), 2196 (2004). 
[37] P. E. Barclay, K. Srinivasan, M. Borselli, and O. Painter, Appl. 

Phys. Lett. 85(1), 4 (2004). 
[38] H. M. Gibbs, Optical testability: controlling light with light 

(Academic Press, Orlando, FL, 1985). 
[39] M. Borselli, K. Srinivasan, P. E. Barclay, and O. Painter, Appl. 

Phys. Lett. 85(17), 3693 (2004). 
[40] P.-T. Lee, J. R. Cao, S.-J. Choi, Z.-J. Wei, J. O'Brien, and P. D. 

Dapkus, Appl. Phys. Lett. 81(18), 3311 (2002). 
[41] A. Cutolo, M. Iodice, P. Spirito, and L. Zeni, J. Lightwave Tech. 

15(3), 505 (1997). 
[42] G. Cocorullo and I. Rendina, IEE Elec. Lett. 28(1), 83 (1992). 
[43] G. P. Agrawal and N. K. Dutta, Semiconductor Lasers (Van 



Nostrand Reinhold, New York, NY, 1993). 
[44] S. G. Johnson, M. Ibanescu, M. A. Skorobogatiy, O. Weis- 
berg, J. D. Joannopoulos, and Y. Fink, Phys. Rev. E 65, 06661 1 
(2002). 

[45] One can show l43l l44tl that the appropriate weighting is 
« 2 (r)£' 2 (r), proportional to the local electric field energy den- 
sity. Note also that the time averaged electromagnetic energy 
stored in the cavity, U, is J (1 /2)e n 2 (r)E 2 (r)dr for a harmonic 
mode in a non-dispersive dielectric. 

[46] We neglect the small (< 1 x 10 16 cm~ 3 ) background free- 
carrier hole density due to the ionized acceptors of the p-type 
Si layer used in this work. 

[47] In using the approximate theory above, in which regions of high 
two-photon absorbed power are correlated with high steady- 
state carrier density, we better approximate the cavity "volume" 
of interest, and consequently the effective free-carrier lifetime 
better represents the average time a free-carrier stays in the re- 
gion of the PC cavity mode. 

[48] This expression neglects the tensor nature of the third order sus- 
ceptibility. A more general expression, which includes the ten- 
sor nature of the Kerr effect, is given in Ref. 1 9]. 

[49] Experimental results| 30] indicate that this Drude model must be 
modified slightly to accurately describe the contribution from 
hole free-carriers in silicon, which scales with A^' 8 . For sim- 
plicity, we ignore this in the following analysis, and note that 
the modification is straightforward, and is included in later nu- 
merical results. 

[50] Heating due to each absorption process will result in a spatial 
temperature profile given by solving the heat equation. In gen- 
eral the different processes have different spatial heating, and 
therefore temperature, profiles. 

[51] Note that the sharp transition edge associated with optical bista- 
bility occurs at the cavity resonance wavelength when scanning 
from blue to red, thus an accurate measure of AX„ can be made, 
even if one of 5A, cannot. 

[52] Note that since we have two data points for each input power, 
one dispersive and one absorptive, there will be an optimum 
x(f,) with non-zero residual error. 



